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Abstract 

A new Lagrangian model without nonlinear scalar self-interactions in the relativistic mean- 
field (RMF) theory is proposed. Introducing terms for scalar- vector interactions (SVI), we have 
developed a RMF Lagrangian model for finite nuclei and nuclear matter. It is shown that by 
inclusion of SVI in the basic RMF Lagrangian, the nonlinear and ci^ terms can be dispensed 
with. The SVI Lagrangian thus obtained provides a good description of ground-state properties 
of nuclei along the stability line as well as far away from it. This Lagrangian model is also able to 
describe experimental data on the breathing-mode giant monopole resonance energies well. 

PACS numbers: 21.30.Fe, 21.10.Dr, 21.60.-ii, 24.10.Cii, 24.10.Jv, 2im.+i 
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The relativistic mean- field (RMF) theory [l|, l2|, l3| of nuclear interaction based upon 
exchange of mesons between nucleons has established itself as a successful approach to de- 
scribing properties of nuclei along the stability line as well as far away from it 0, ^, Q]. 
The Dirac-Lorentz structure of nucleons provides a built-in spin-orbit interaction with its 
advantage over the non-relativistic Skyrme theory in describing properties such as anoma- 
lous isotope shifts in Pb nuclei j?], which depend upon shell structural effects. An isospin 
dependence of the spin-orbit interaction or rather a lack of it is responsible for the anoma- 
lous behaviour of the isotope shifts Sj . The idea of pseudospin symmetry in nuclei has been 
attributed to the relativistic (Dirac) nature of nucleons 9|. 

n 

With the advent of the Walecka model for nuclei and nuclear matter, it was realized 

that the linear RMF Lagrangian with its large value of nuclear incompressibility was unable 

to describe properties of finite nuclei properly. A lack of proper ingredients for a suitable 

description of surface properties was cited as a main reason for this drawback. In order 

□ 

to remedy this problem, Boguta and Bodmer pj| introduced nonlinear scalar self-coupling 
terms of the form + a^. Consequently, the RMF Lagrangian has proved to be successful 
and the nonlinear scalar terms have thus become indispensable for an adequate treatment 
of finite nuclei and nuclear matter saturation. The nonlinear cr-field seems to provide an 
essential density dependence of nuclear force in a finite nucleus. With the inclusion of the 
nonlinear cx-field, the theory also becomes renormalizable. 

One of the first successful nuclear forces within this model is NL-SH jsj. Within this 
model an improved set NL3 has been obtained for finite nuclei Gj. However, as is well 
known, the model with + gives an equation of state (EOS) of nuclear matter that is 
very stiff and is consequently untenable for the spectrum of observed neutron star masses. 

The quartic vector coupling of the form cj^ was introduced in the RMF Lagrangian. 
An appropriate description of finite nuclei was obtained with the vector self-coupling of uj 
meson [121, 13] • The inclusion of the cj^ coupling has also helped to improve the shell effects 
along the stability line |l3|. A desired softe ning of the EOS of nuclear matter due to vector 
self-coupling of uo meson was shown in ref. 14 1. 

The density dependence of the nuclear interaction in the RMF theory remains an open 
problem. The Walecka model offers an appropriate avenue to construct a theory that should 
be suitable for describing various aspects of finite nuclei all along the periodic table (an 
ambitious goal) as well as properties of nuclear matter concomitantly. Point-coupling models 



2 



have been introduced 
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IGj to describe finite nuclei. Attempts have been made to broaden 



the basis of the RMF Lagrangian by including terms of higher orders in the scalar and 
vector fields with inclusion of interaction terms amongst various mesonic degrees of freedom 
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have been introduced 



with a view to modify density dependence of the nuclear interaction in an explicit form 
with a good degree of success. This requires inclusion of additional parameters to model 
density-dependence of meson couplings. Notwithstanding the above, the RMF theory serves 
as an ideal platform for an effective field theoretical approach for many-body problems of 
nuclei with sufficient space for innovation. 

An upsurge in experimental data especially in the domain of extreme regions of the 
periodic table provides an incentive to devise new and improved approaches and models 
to be able to describe the same. Savushkin et al. have incorporated various meson- 
meson interactions in their approach especially those between a and u meson in addition 
to nonlinear couplings of both these mesonic fields. This problem has been approached 13] 
from a more general point of view by taking expansion in and interactions amongst various 
mesonic fields. This approach has led to an improvement for finite nuclei and nuclear matter 
with a larger number of parameters required. 

In this work, we have sought to explore the possibility of dispensing with the nonlinear 
scalar self-couplings which have so far remained essential for finite nuclei. We ask ourselves: 
whether it is possible to mock the scalar self-couplings and their inherent density dependence 
in nuclei by employing meson-meson interactions especially between a and u mesons instead? 
Keeping the issue of renormalizability in abeyance, we have added couplings between a and 
u mesons of the form au'^ + a^cu^ to the basic (linear) RMF Lagrangian based upon exchange 
of a, uj and p mesons. Properties of nuclear matter for the scalar-vector interaction (SVI) 
of this form were explored in ref. 2^. Recently, the Lagrangian model SIG-OM with the 
inclusion of the coupling of the form a'^uo'^ whilst retaining the scalar self-couplings + o"^ 
has been developed |26j] . In the present work, we have narrowed down the space by excluding 
the self-couplings at the expense of the scalar-vector meson-meson couplings. 

The basic RMF Lagrangian density that describes nucleons as Dirac spinors interacting 
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with the meson fields is given by 



-\r,.R''' + \mlp,f?^ - \f,^F^\ (1) 

where M^v is the bare nucleon mass and ip is its Dirac spinor. Nucleons interact with a, uo, 
and p mesons, with couphng constants being g^, g^ and (jfp, respectively. The photonic field 
is represented by the electromagnetic vector A'^. The effective Lagrangian for finite nuclei 
that is used commonly is given by 

-Ceff = -^^0 ~ f^NL- (2) 

The nonlinear cr-meson self-couplings which have so far been an integral part of the RMF 
Lagrangian are given by 

U^L = \92(y^ + \gz<T^. (3) 

The parameters g2 and (73 are the nonlinear couplings of a-meson in the conventional + 
a'^ model 1^. Here, we put g2 = gs = 0, thus eliminating the self-couplings [/nl of a meson. 
Instead, we introduce the meson-meson interaction terms of the form 

where (74 and g^ represent the respective coupling constants for meson-meson interactions 
between a and u mesons. The effective Lagrangian in our case is then 

The corresponding Klein-Gordon equations can be written as 

(-A + mf)(T = -g^^Pij 
{-A + m*J)uj^ = g^il)^yilJ 
(-A + mJ)p; = gpip-fuTifj 

-AA, = ^67/^(1 + r3)7.^, (6) 
where the effective meson masses m* and m* can be obtained as 

*2 2 2 lln \ 2 

= m^- g^ujo / {2a) - g^u^ 

= ml + g^a + g^a'^. (7) 



These equations represent an implicit density dependence of a and uj meson masses and 
effectively that of the nuclear interaction therein. 

The parameters of the new Lagrangian model SVI are obtained by a multi-dimensional 
search in the parameter space by fitting experimental binding energies and charge radii of 
a set of a few nuclei (cf. |5|] for a detailed procedure). The nuclei included are ^^O, ^°Ca, 
^°Zr, ^^^Sn, ^^^Sn and ^°^Pb. The isotopes ^^^Sn and ^^"^Sn are included in order to span the 
broad range of isospin. No conditions have been put on nuclear matter properties and thus 
the parameters are allowed to vary freely without any bias to the nuclear matter properties. 
The uo and p meson masses have been fixed at their empirical values. 

The parameters of the Lagrangian obtained as a result of a free variation in the multi- 
dimensional space are shown in Table I. We have obtained two parameter sets SVI-1 and 
SVI-2 which are deemed as appropriate for ground-state binding energies and charge radii 
of nuclei. The parameter of the forces NL-SH and NL3 with the nonlinear scalar couplings 
are also shown for comparison. 

The nuclear matter properties of SVI-1 and SVI-2 are shown in the lower section of Table 

I. The sets SVI-1 and SVI-2 are close to each other in the nuclear matter properties with a 
slight difference in the incompressibility with K = 264 MeV for SVI-1 and K = 272 MeV 
for SVI-2. There is only a minor difference in the effective mass m*. The m* values for SVI 
interactions are clearly higher than those of the Lagrangian sets NL-SH and NL3 with the 
scalar self-interactions. 

The saturation density for both SVI-1 and SVI-2 is slightly higher than that of NL-SH 
and NL3. One notable difference between the two SVI sets is the difference in the asymmetry 
energy J (or 04). One can note that even in an interaction that is different from the scalar 
self-interactions, it is not possible to bring down the asymmetry energy in the acceptable 
range of 30-33 MeV. 

The binding energy of key spherical nuclei along with a few representative ones as obtained 
with SVI-1 and SVI-2 is shown in Table II. For a comparison, we also show the results due 
to NL-SH and NL3. With the exception of the very light nucleus of ^^O, both SVI-1 and 
SVI-2 show an excellent agreement with the experimental binding energies over a large range 
of mass. This is reflected by the smaller value of the rms deviation 6 of SVI-1 and SVI-2 
values from the experimental data vis-a-vis NL3 and NL-SH shown at the bottom of Table 

II. A marked improvement with SVI interactions is in the binding energy of doubly magic 
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TABLE I: The parameters and nuclear matter (NM) properties of the scalar-vector Lagrangians 
SVI-1 and SVI-2 without nonlinear scalar self-couplings. The sets NL-SH and NL3 with the scalar 
self-couplings are also shown for comparison. 



Parameters 


SVI-1 


SVI-2 


NL-SH 


NL3 


M (MeV) 


939.0 


9.39.0 


939.0 


939.0 


rua- (MeV) 


524.527 


524.024 


526.0592 


508.194 


(MeV) 


783.0 


783.0 


783.0 


782.501 


mp (MeV) 


763.0 


763.0 


763.0 


763.0 




9.6762 


9.641 


10.4436 


10.217 




11.6028 


11.565 


12.9451 


12.8675 


9p 


4.464 


4.492 


4.3828 


4.4744 


§2 (fm-^) 








-6.9099 


-10.432 










-15.8337 


-28.885 


54 (fm^i) 


17.1537 


16.962 


0.0 


0.0 


55 


33.8565 


32.819 


0.0 


0.0 


NM properties 










Po (fm-3) 


0.149 


0.149 


0.146 


0.148 


ay (MeV) 


-16.30 


-16.31 


-16.33 


-16.24 


K (MeV) 


263.9 


271.5 


354.9 


271.6 


m* 


0.616 


0.621 


0.597 


0.595 


J (MeV) 


37.6 


37.0 


37.0 


37.4 



nuclei ^°°Sn, ^^^Sn and ^°^Pb over those of NL-SH and NL3. This may have consequences 
on the shell effects in nuclei especially in the vicinity of the r-process path and drip lines. 

The charge radii of nuclei obtained with SVI-1 and SVI-2 are shown in Table III. These 
are compared with the values obtained with NL-SH and NL3. The SVI interactions describe 



the experimental data [27[ on nuclei well. An improvement on the charge radii of Pb isotopes 
over those of NL3 can be seen. 

We have calculated the ground-state properties of the isotopic chains of Sn and Pb. 
Especially, the chain of Sn isotopes offers experimental binding energies over the whole 
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TABLE II: The binding energy (in MeV) of nuclei calculated with SVI-1 and SVI-2 and compared 
with NL3 and NL-SH. The empirical values (exp.) are shown in the last column. The rms deviation 
(6) of theoretical values from the experimental data is shown at the bottom of the table. 



Nucleus 


SVI-1 


SVI-2 


NL-SH 


NL3 


exp. 


160 


-129.7 


-129.7 


-128.4 


-128.8 


-127.6 


40Ca 


-343.2 


-343.2 


-340.1 


-342.0 


-342.1 


48Ca 


-415.3 


-415.0 


-415.1 


-415.2 


-416.0 


902J. 


-783.0 


-783.0 


-782.9 


-782.6 


-783.9 


looSn 


-827.3 


-827.1 


-830.6 


-829.2 


-824.5 


ii^Sn 


-988.3 


-988.4 


-987.9 


-987.7 


-988.7 


i24Sn 


-1049.7 


-1049.6 


-1050.1 


-1050.2 


-1050.0 


i32Sn 


-1103.8 


-1103.3 


-1105.9 


-1105.4 


-1102.9 


202pb 


-1591.2 


-1591.8 


-1595.8 


-1592.6 


-1592.2 


208pb 


-1637.0 


-1637.3 


-1640.4 


-1639.5 


-1636.7 




-1661.7 


-1662.3 


-1664.3 


-1661.6 


-1663.3 


6 


1.33 


1.20 


2.70 


2.00 





TABLE III: The rms charge radius (in fm) obtained with SVI-1 and SVI-2. The values for NL-SH 
and NL3 are also shown. The rms deviation (6) of theoretical values from the experimental data 
is shown at the bottom of the table. 



Nucleus 


SVI-1 


SVI-2 


NL-SH 


NL3 


exp. 


160 


2.698 


2.700 


2.699 


2.728 


2.730 


^"Ca 


3.438 


3.442 


3.452 


3.470 


3.450 


4«Ca 


3.451 


3.457 


3.452 


3.470 


3.450 


90Zr 


4.277 


4.285 


4.289 


4.287 


4.258 


ii^Sn 


4.593 


4.601 


4.599 


4.599 


4.626 


i24Sn 


4.644 


4.652 


4.651 


4.661 


4.673 


208pb 


5.501 


5.508 


5.509 


5.523 


5.503 


214pb 


5.562 


5.568 


5.562 


5.581 


5.558 


S 


0.021 


0.019 


0.021 


0.020 
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range from the doubly magic nucleus ^^'^Sn to the doubly magic ^^^Sn, thus encompassing 
the space between the two magic numbers = 50 and = 82. The difference AE of 
binding energy of nuclei obtained with RMF+BCS calculations from the experimental value 
is shown for Sn and Pb isotopes in Fig. 1. For the BCS pairing, neutron pairing gaps have 
been obtained from the experimental masses of neighbouring nuclei. 

The problem of the arches and a predominance of shell energy at the magic numbers is 
well-known. It pervades both the microscopic theories as well as macroscopic-microscopic 
mass formulae. Viewing the results for Sn isotopes in Fig. 1(a), one can notice unambiguous 
arches at the two magic numbers especially with the force NL3. With exception of the region 
near ^''^Sn, both SVI-1 and SVI-2 describe the data well for nuclei including those near ^^^Sn 
(A^ = 82). For nuclei in the vicinity of ^^''Sn (A^ = 50) SVI-1 and SVI-2 show a significant 
improvement over the results of NL3. The rms deviation from the experimental data is 1.04 
MeV and 1.07 MeV for SVI-1 and SVI-2, respectively. In comparison, it is 1.83 MeV for 
NL3. This indicates a significant improvement in the binding energies for Sn isotopes with 
SVI-1 and SVI-2, especially near the magic numbers A^ = 50 and A^ = 82. The arch-like 
behaviour with SVI-1 and SVI-2 is reduced considerably. 

This pattern is also visible for the isotopic chain of Pb in Fig. 1(b). Both SVI-1 and SVI-2 
exhibit a significant improvement in the binding energies over NL3. The rms deviation of 
the theoretical values with SVI-1 and SVI-2 is 0.88 MeV and 0.69 MeV, respectively. This is 
much smaller than the corresponding value of 1.82 MeV with NL3 for the Pb isotopes. Thus, 
SVI interactions provide a better description of the binding energies of the Pb isotopes. In 
comparison, NL3 values overestimate the data near the magic number and gives a well- 
formed arch about the magic number. With NL3 divergences of the binding energies near 
the magic number are displayed strongly as has been observed also for the Sn isotopes above. 
It is a matter of further investigation as to what ingredients in the RMF theory would lead 
to divergences or a lack thereof at shell closures. 

The charge radii of Pb isotopes and the anomalous kink in charge radii represent a 
characteristic feature related to shell structure of nuclei. It was shown that the RMF theory 
with NL-SH reproduced the anomalous kink successfully [7]. This feature has since been 
demonstrated by all the Lagrangian models in the RMF theory. In order to discern the 
behaviour of SVI in this respect, we show in Fig. 2 the charge radii of Pb isotopes calculated 
with SVI-1 and SVI-2. The force SVI-1 reproduces the experimental data 271] over the range 



8 



> 

CD 



LU 

< 



2 

-2 
-4 
-6 



> 

0) 



LU 

< 



I ' ' ' I ' ^ 

Sn 



~l 1 T" 



~l 1 T" 



(a)J 
o 



» • SVI-1 

0---0 SVI-2 
■ ■ NL3 



_i I i_ 



_l I L. 



100 



108 



116 



124 



132 




216 



FIG. 1: Binding energies for (a) Sn and (b) Pb isotopes with SVI-1 and SVI-2 and compared with 
the experimental data. The results obtained with NL3 are also shown. 

of Pb isotopes very well. In comparison, SVI-2 values overestimate the data only slightly. 
The set NL3, on the other hand, overestimates the charge radii of all the Pb isotopes 
significantly. 

In order to test the applicability of the new model for nuclei away from the line of (3- 
stability, we have performed axially deformed RMF calculations for a number of nuclei. 
The nuclei encompass representative cases from ^^Si to ^^^Pt, which includes nuclei from 
medium masses through the rare-earth region to higher masses. The results of calculations 
with SVI-1 and SVI-2 are shown in Table IV. The experimental data on the binding energies 
are taken from the recent high-precision mass measurements on Si 281], Sr 291] and Mo 30]. 
The binding energy of other nuclei has been taken from the 2003 compilation of atomic 
masses (3l|. 



Both SVI-1 and SVI-2 provide an excellent description of the experimental binding en- 
ergies of nuclei over a large range of atomic mass. The difference between the predictions 
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FIG. 2: Charge radii of Pb isotopes obtained with SVI-1, SVI-2 and NL3 and compared with the 
experimental data. 

of the two sets are small. For a few cases SVI-1 provides a better description whereas for a 
few others the SVI-2 does better. The rms deviation of the binding energies for both the 
sets amounts to ~ 0.62 MeV. 

The original Walecka model (linear) [l] with the nucleon- meson couplings of a and lo 
mesons has been instructive for achieving saturation of nuclear matter. With the inclusion 
of nonlinear scalar self-couplings, the saturation is achieved at nuclear matter properties 
viz., the compressibility within the acceptable range (cf. nuclear matter properties in Table 
I with e.g. NL-SH and NL3). Here the incompressibility of nuclear matter K denotes the 
cardinal point on the saturation curve. The incompressibility K plays an important role 
in determining the breathing-mode giant monopole resonance (GMR) energies. It is thus 
important that an acceptable Lagrangian model should be able to describe the GMR data. 

For a comparative analysis of the Lagrangian models involved, we have carried out con- 
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TABLE IV: The binding energy (in MeV) and quadrupole deformation (32 (in parentheses) of 
nuclei away from the stabihty Une calculated with SVI-1 and SVI-2. The experimental values 
(exp.) where available are shown for comparison. 



Nucleus 


SVI-1 


SVI-2 


exp. 


36 Si 


-292.8 (-0.03) 


-292.7 (0.03) 


-292.0 


38 Si 


-300 6 fO 28) 


-300 2 (0 28) 


-299.9 


40si 


-306.9 (0.35) 


-306.4 (0.35) 


-306.5 


40s 


-333 2 fO 25) 


-332 8 (0 25) 


-333.2 


76^1 


-632.9 (0.0) 


-632.3 (0.0) 


-633.1 


86sr 


-748.5 (0.0) 


-748.5 (0.0) 


-748.9 


88si. 


-768.1 (0.0) 


-768.0 (0.0) 


-768.5 


hoMo 


_918.6 (-0.23) 


-918.5 (-0.23) 


-919.5 


i20Xe 


-1008.1 (0.28) 


-1008.4 (0.27) 


-1007.7 (0.22) 


i40Xe 


-1161.3 (0.10) 


-1161.2 (0.10) 


-1160.7 (0.11) 


i54Sm 


-1267.4 (0.31) 


-1267.5 (0.30) 


-1266.9 (0.34) 


164Dy 


-1337.9 (0.35) 


-1338.1 (0.35) 


-1338.0 (0.35) 


168 


-1365.3 (0.34) 


-1365.5 (0.34) 


-1365.8 (0.34) 


174Yb 


-1406.7 (0.31) 


-1407.0 (0.31) 


-1406.6 (0.33) 


190^ 


-1510.5 (0.19) 


-1510.7 (0.20) 


-1509.9 


196pt 


-1552.2 (0.12) 


-1552.5 (0.12) 


-1553.6 (0.13) 



strained Generator Coordinate Method (GCM) calculations [321] for the isoscalar GMR mode 
for a set of nuclei. The nuclei included are ^°Zr, ^^°Sn and ^°^Pb. The results of GCM cal- 
culations are shown in Table V. 

The set SVI-1 with K = 264 MeV underestimates the experimental values for ^°Zr and 
^°^Pb by more than ~ 0.5 MeV, whereas for ^^°Sn, the disagreement with the datum is 
nominal. In comparison, SVI-2 with K = 272 MeV gives GMR energies for ^°Zr and ^^''Sn, 
which are very close to the experimental data. For ^°^Pb, its values is slightly smaller than 
the experimental one. Comparatively, NL3 with K = 272 MeV gives values which are 
systematically smaller than those of SVI-2. This difference in predictions of a nonlinear 
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TABLE V: The breathing mode GMR energies obtained with constrained GCM calculations using 
SVI-1 and SVI-2. The experimental data 33, 3^ are also shown. 



Nucleus SVI-1 SVI-2 NL3 exp. 

90Zr 17.2 17.5 16.9 17.81 ± 0.30 

i20Sn 15.2 15.4 15.0 15.52 ±0.15 

208pb 13.3 13.5 13.0 13.96 ±0.28 



scalar Lagrangian from those of a scalar-vector Lagrangian points to some subtle differences 
in finite effects, possibly surface, of the two models. A comparison between the two sets of 
SVI shows that SVI-2 is commensurate with the ground-state properties of finite nuclei as 
well as with the GMR energies. 

In conclusion, we have constructed a Lagrangian model without nonlinear scalar self- 
couplings in the RMF theory. Incorporating scalar-vector interaction terms in the linear 
RMF Lagrangian, the Lagrangian model SVI has been developed. We have obtained pa- 
rameter sets SVI-1 and SVI-2 in the framework of the new model. It is shown that both 
SVI-1 and SVI-2 provide a good description of the ground-state properties of nuclei along 
the stability line as well as for nuclei far away from it. With an incompressibility of nuclear 
matter K ~ 272 MeV, SVI-2 is able to reproduce the breathing-mode GMR energies on key 
nuclei well. Thus, the model SVI without nonlinear interactions becomes viable for finite 
nuclei and nuclear matter. 

I thank Prof. Lev Savushkin for useful discussions. This work is supported by the 
Research Administration Project No. SP07/03 of Kuwait University. 
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